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Abstract. I review the constraints imposed by the observed extragalactic back- 
ground light (EBL) on the history of the stellar birthrate in galaxies. At faint 
magnitudes, the logarithmic slope of the galaxy counts is flatter than 0.4 in all 
seven UBVIJHK optical bandpasses of the Hubble Deep Fteld-South imaging sur- 
vey. The integration of the number counts provides a lower limit to the surface 
brightness of the optical extragalactic sky of 15 nWm -2 sr -1 , comparable to the 
intensity of the far-IR background from COBE data. If the initial mass function has 
a Salpeter slope with a lower mass cutoff consistent with observations of M subdwarf 
disk stars, a lower limit of Q, > 0.005 J50 (at Hubble constant 50 kms -1 Mpc -1 ) is 
derived for the visible (processed gas + stars) mass density needed to generate an 
extragalactic background light (EBL) at a level of 50 /so nWm -2 sr _1 . The current 
'best-guess' estimate to is 0.012 /50, about 16% of the nucleosynthetic baryon 
density. The contribution of quasar activity to the observed EBL is unlikely to 
exceed 20%. 



1 Introduction 

Recent progress in our understanding of faint galaxy data made possible by 
the combination of Hubble Space Telescope (HST) deep imaging and ground- 
based spectroscopy has dramatically increased our knowledge of the evolution 
of the stellar birthrate in optically-selected galaxies from the present-epoch 
up to z as 4 p2| , |3j|] , f2]]| . The explosion in the quantity of information 
available on the high-redshift universe at optical wavelengths has been com- 
plemented by the measurement of the far-IR/sub-mm background by DIRBE 
and FIRAS onboard the COBE satellite ^ , || , §7j] , by the detection of dis- 
tant ultraluminous sub- mm sources with the SCUBA camera p8| , Q , and by 
theoretical progress made in understanding how intergalactic gas follows the 
dynamics dictated by dark matter halos until radiative, hydrodynamic, and 
star formation processes take over [^j, |p0| , The IR data have revealed 
the 'optically-hidden' side of galaxy formation, and shown that a significant 
fraction of the energy released by stellar nucleosynthesis is re-emitted as ther- 
mal radiation by dust , |l8| . The underlying goal of all these efforts is to 
understand the growth of cosmic structures, the internal properties of galax- 
ies and their evolution, and ultimately to map the star formation history of 
the universe from the end of the cosmic 'dark age' to the present epoch. 

In this talk I will focus on the galaxy number-apparent magnitude relation 
and its first moment, the integrated galaxy contribution to the extragalac- 
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tic background light (EBL). The logarithmic slope of the differential galaxy 
counts is a remarkably simple cosmological probe of the history of stellar birth 
in galaxies, as it must drop below 0.4 to yield a finite value for the EBL. The 
recently released Hubble Deep Field-South (HDF-S) images, together with 
other existing HST and ground-based observations, provide a unique dataset 
to estimate the spectrum and amplitude of the optical EBL from discrete 
sources. Together with the far-IR/sub-mm background, the optical EBL is 
an indicator of the total luminosity of the universe, as the cumulative emis- 
sion from young and evolved galactic systems, as well as from active galactic 
nuclei (AGNs), is recorded in this background. As such it provides, for a given 
initial mass function, a quantitative estimate of the baryonic mass that has 
been processed by stars throughout cosmic history. 

Unless otherwise stated, an Einstein-de Sitter (EdS) cosmology {Qm = L 
Qa = 0) with H = 100 /i kins" 1 Mpc" 1 will be adopted in this talk. All 
magnitudes will be given in the AB system. The work presented here has 
been done in collaboration with L. Pozzetti. 



2 Galaxy counts 

The HDF-S dataset includes deep near-IR NICMOS images and the deepest 
observation ever made with the STIS 50CCD filtcrless imaging mode. The 
galaxy sample used here was extracted from version 1 of the HDF-S catalog 



on ftp:/ /archive , stsci.edu/pub/hdf_south/versionl/. At near-IR wavelengths 
(in the F110W, F160W, and F222M bandpasses, corresponding to the J, H, 
and K filters), it consists of 425 objects detected in the J + H image, over a 
field of 50" x 50". The 50CCD (corresponding roughly to a V + I filter) STIS 
catalog consists of 674 objects detected again over a field of the same size. 

Figure 1 shows the HDF-N and -S galaxy counts compiled directly from 
the catalogs, as a function of AB isophotal magnitudes in the UBVIJHK 
bandpasses for all galaxies with signal-to-noise ratio S/N > 3 within the 
band. No correction for detection completeness have been made. A compi- 
lation of existing HST and ground-based data is also shown [p5| . All 
magnitudes have been corrected to the AB system, while the second order 
colour corrections for the differences in the filter effective wavelengths have 
not been applied to the ground-based data (for the typical colours of galax- 
ies in the HDF these corrections are less than 0.1 mag). The HDF optical 
counts agree well with previous surveys, to within 20% in the magnitude 
range 22 < m^B < 26. One should note, however, that different algorithms 
used for 'growing' the photometry beyond the outer isophotes of galaxies 
can significantly change the magnitude of faint galaxies. According to Q, 
roughly 50% of the flux from resolved galaxies with V > 23 mag lie outside 
the standard-sized apertures used by photometric packages. An extragalactic 
sky pedestal created by the overlapping wings of resolved galaxies may also 
contribute significantly to the sky level, and would be undetectable except 
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Fig. 1. Left: Differential UBVIJHK galaxy counts as a function of AB magni- 
tudes. The sources of the data points are given in the text. Note the decrease of 
the logarithmic slope dlogN/dm at faint magnitudes. Right: Extragalactic back- 
ground light per magnitude bin, i v = 10~°' 4(mAB+48 6) 7V(m), as a function of U 
{filled circles), B {open circles), V {filled pentagons), I {open squares), J {filled 
triangles), H {open triangles), and K {filled squares) magnitudes. For clarity, the 
BVIJHK measurements have been multiplied by a factor of 2, 6, 15, 50, 150, and 
600, respectively p^ |. 

by absolute surface photometry ||. Also, at faint magnitude levels, distant 
objects which are brighter than the nominal depth of the catalog may be 
missed due to the (1 + z) 4 dimming factor. All these systematic errors are 
inherent in HST faint-galaxy photometry; as a result, our estimate of the 
integrated flux from resolved galaxies will typically be too low, and must be 
strictly considered as a lower limit. 



3 The extragalactic background light 

The contribution of known galaxies to the optical EBL can be calculated 
directly by integrating the emitted flux times the differential number counts 
down to the detection threshold. I have used the compilation of ground- 
based, HST, and HDF data shown in Figure 1 to compute the integrated flux 
at 0.35 < A < 2.2 /im. 

In all seven bands, the slope of the differential number-magnitude relation 
is flatter than 0.4 above m^B ~ 20 (25) at near-IR (optical) wavelengths, and 
this flattening appears to be more pronounced at the shorter wavelengths. 
The leveling off of the counts is clearly seen in Figure 1, where the func- 
tion i v = io~°- 4 ( mAB + 48 - 6 ) N{m) is plotted against apparent magnitude in all 
bands. While counts having a logarithmic slope dlogN/dniAB = ct > 0.40 
continue to add to the EBL at the faintest magnitudes, it appears that the 
HDF survey has achieved the sensitivity to capture the bulk of the near- 
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ultraviolet, optical, and near-IR extragalactic light from discrete sources. The 
flattening at faint apparent magnitudes cannot be due to the reddening of 
distant sources as their Lyman break gets redshifted into the blue passband, 
since the fraction of Lyman-break galaxies at (say) B ks 25 is only of order 
10% [^2|. Moreover, an absorption-induced loss of sources cannot explain the 
similar change of slope in the V, I, J, H, and K bands. While this suggests 
that the surface density of optically luminous galaxies is leveling off beyond 
z ~ 1.5, one should worry about the possibility of a significant amount of 
light being missed at faint magnitudes. 



log M/Am) 



Fig. 2. Spectrum of the optical extragalactic background light from resolved sources 
as derived from a compilation of ground-based and space-based galaxy counts in 
the UBVIJHK bands {filled dots), together with the FIRAS 125-5000 fan {solid 
and dashed lines) and DIRBE 140 and 240 /im {filled squares) detections ]is[ |, p|. 
The empty squares show the DIRBE points after correction for WIM dust emission 
@. Also plotted {filled triangle) is a FOCA-UV point at 2000 A (l[, and a tentative 
detection at 3.5 fim {empty dot) from CCRE/DIRBE observations Q. The empty 
pentagons at 3000, 5500, and 8000 A are recent detections from absolute photometry 
Q. Upper limits are from jli), the lower limit from Jl(J. Note that the values 
obtained by integrating the brightness of resolved galaxies are strict lower limits to 
the EBL intensity. 



The spectrum of the optical EBL is shown in Figure 2, together with 
the recent results from COBE. The value derived by integrating the galaxy 
counts down to very faint magnitude levels (because of the flattening of 
the number-magnitude relation most of the contribution to the optical EBL 
comes from relatively bright galaxies) implies a lower limit to the EBL in- 
tensity in the 0.2-2.2 /jm interval of J op t ~ 15 nWm~ 2 sr _1 . Including the 
tentative detection at 3.5 /zm by |?J would boost I pt to w 19 nWm~ 2 sr _1 . 
Recent estimates of the optical EBL at 3000, 5500, and 8000 A from abso- 
lute surface photometry by || lie between a factor of two to three higher 
than the integrated light from galaxy counts. Applying this correction fac- 
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tor to the range 3000-8000 A gives an optical EBL intensity in excess of 
25nWm- 2 sr 1 in the interval 0.2-3.5 fxm. The CORE/FIRAS |1| mea- 
surements yield /fir ~ 14 nWm _2 sr _1 in the 125-2000 /im range. When 
combined with the DIRBE |l9), |27), points at 140 and 240 /im, one gets 
a far-IR background intensity of 7 F ir(140 - 2000 jum) « 20 nWm^sr 1 . 
The residual emission in the 3.5 to 140 fim region is poorly known, but it is 
likely to exceed 10 nWm~ 2 sr -1 ||. A 'best-guess' estimate to the total EBL 
intensity observed today is then 

/ EB L = 55±20nWm" 2 sr 1 . (1) 

In the rest of my talk, I will adopt a reference value for the background light 
associated with star formation activity over the entire history of the universe 
of 50/ 5 onWm- 2 sr- 1 . 



4 EBL from quasar activity 

A direct estimate of the contribution of quasars to the EBL depends on the 
poorly known bolometric correction and the possible existence of a distant 
population of dusty AGNs with strong intrinsic absorption, as invoked in 
many models for the X-ray background. These Type II QSOs, while unde- 
tected at optical wavelengths, could contribute significantly to the far-IR 
background. It is in principle possible to bypass some of the above uncer- 
tainties by weighing the local mass density of black holes remnants j29j. 
Recent dynamical evidence indicates that supermassive black holes reside at 
the center of most nearby galaxies. The available data (see Fig. 3) show a 
correlation (but with a large scatter )between bulge and black hole mass, 
with Mbh ~ 0.006 M sp h as a best-fit The mass density in old spheroidal 
populations today is estimated to be £2 S phh = 0.0018io]ooo85 iW i implying a 
mean mass density of quasar remnants today 

Pbh = 3 ± 2 x 10 6 h Mq Mpc -3 . (2) 

Since the observed energy density from all quasars is equal to the emitted 
energy divided by the average quasar redshift, the total contribution to the 
EBL from accretion onto black holes is 

/bh = Y~TT^T ~ 4 ± 2.5 nWm~ 2 sr _1 770.0577-7— r (3) 

47T (1 + Z) (1+ 2 ) 

(h = 0.5), where 170.05 is the efficiency of accreted mass to radiation conversion 
(in units of 5%). Therefore, unless dust-obscured accretion onto supermassive 
black holes is a very efficient process (77 3> 0.05), a population of quasars 
peaking at z ~ 1.5 — 2 is expected to make a contribution to the brightness 
of the night sky not exceeding 10-20% @, (2l|. 
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Fig. 3. Black hole mass distribution against the bulge luminosity of their host 
galaxies P3| . Arrows indicate upper limits on Mbh- The symbols denote different 
galaxy types: empty circles (E), filled squares (SO), filled circles (Sab), and empty 
squares (Sbc-Scd). 



5 The stellar mass density today 



With the help of some simple stellar population synthesis tools we can now 
set a lower limit to the total stellar mass density that produced the observed 
sky brightness and constrain the cosmic history of star birth in galaxies. One 
of the most serious uncertainties in this calculation has always been the lower 
cutoff (usually treated as a free parameter) of the initial mass function (IMF). 
Observations of M subdwarfs stars with the HST have recently shed some 
light on this issue, showing that the IMF in the Galactic disk can be repre- 
sented analytically over the mass range 0.1 < m < 1.6 (here m is in solar 
units) by log</>(m) = const - 2.33 logm - 1.82(logm) 2 (@, hereafter GBF). 
For m > 1 this mass distribution agrees well with a Salpeter function. Obser- 
vations of normal Galactic star-forming regions also show some convergence 
in the basic form of the IMF at intermediate and high masses, a power-law 
slope that is consistent with the Salpeter value |l3| ]. In the following I will 
use a 'universal' IMF with the GBF form for m < 1, matched to a Salpeter 
slope for to > 1; the mass integral of this function is 0.6 times that obtained 
extrapolating a Salpeter function down to 0.1 M Q .[] 

As shown in Figure 4, the bolometric luminosity as a function of age 
r of a simple stellar population (a single generation of coeval, chemically 
homogeneous stars having total mass M, solar metallicity, and the above 



1 The bolometric light contributed by stars less massive than 1 Mq is quite small 
for a 'typical' IMF. The use of the GBF mass function at low masses instead 
of Salpeter leaves then the total radiated luminosity of the stellar population 
virtually unaffected. 
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Fig. 4. Left: Synthetic J4| bolometric luminosity versus age of a simple stellar popu- 
lation having total mass M = 1 Mq, metallicity Z = Z@ (solid line) and Z = 0.2 Z@ 
(dotted line), and a GBF+Salpeter IMF (see text for details). Right: EBL observed 
at Earth from the instantaneous formation at redshift 2f of a stellar population 
having the same IMF, solar metallicity, and mass density Q,h 2 = 0.00175, 0.00125, 
and 0.00075, as a function of zf- Solid curves: EdS universe with h — 0.5. Dashed 
curves: yl-dominated universe with Om = 0.3, Qa = 0.7, and h = 0.65. 



IMF) can be well approximated by 

( 1200L Q ^- t< 2.6 Myr; 

L(r) = 0.7 L Q ^ (^j ^ 2.6 < r < 100 Myr; (4 ) 

l 2-0^01^ (tT^f)" ' 8 t>H)0 Myr 

(cf |^). Over a timescale of 13 Gyr (the age of the universe for an EdS 
cosmology with h = 0.5), about 1.3 MeV per stellar baryon are radiated away. 
This number depends only weakly on the assumed metallicity of stars. In a 
stellar system with arbitrary star formation rate per comoving cosmological 
volume, p*, and formation epoch tp, the comoving bolometric emissivity at 
time t is given by the convolution integral 

pbo\{t) = / L(r)p*(t ~ r)dr. (5) 
Jo 

The total background light observed at Earth (t = f#) is 

w _ ± f ftaM*. (6) 

4tt J tF 1 + z 

where the factor (1 + z) at the denominator is lost to cosmic expansion when 
converting from observed to radiated (comoving) luminosity density. To set 
a lower limit to the present-day mass density, .!?*, of processed gas + stars 
(in units of the critical density p cr i t = 2.77 x 10 11 h 2 M Q Mpc~ 3 ), consider 
now a scenario where all stars are formed instantaneously at redshift zp. 
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The background light that would be observed at Earth from such an event is 
shown in Figure 4 as a function of z F for Q^h 2 = 0.00075,0.00125,0.00175 
(corresponding to 9, 6.5, and 4 percent of the nucleosynthetic baryon den- 
sity, Q\yh 2 = 0.0193 ± 0.0014 ||), and two different cosmologies. A couple of 
points are worth noting here: (1) the time evolution of the luminosity radi- 
ated by a simple stellar population makes the dependence of the observed 
EBL from zp much shallower than the (1 + zp) lost to cosmic expansion, 
as the energy output from stars is spread over their respective lifetimes; and 
(2) in order to generate an EBL at a level of 50 /50 nWm _2 sr _1 , one re- 
quires f2*h 2 > O.OOI25/50 (for an EdS universe with h = 0.5), hence a mean 
mass-to-blue light ratio today of (M/Lb)* > 3.5/50 for a present-day blue 
luminosity density of ps = 2.0 x 10 8 /iL Q Mpc~ 3 [[ll]. As shown in Figure 4, 
the dependence of these estimates on the cosmological model is rather weak. 
With the adopted IMF, about 30% of this mass will be returned to the inter- 
stellar medium in 10 s yr after intermediate-mass stars eject their envelopes 
and massive stars explode as supernovae. This 'return fraction', R, becomes 
50% after about 10 Gyr. 

A visible mass density at the level of the above lower limit, 0.00125ft- -2 /50, 
while able to explain the measured sky brightness, requires all the stars that 
give origin to the observed light to have formed at zp k, 0.2, and is, as such, 
rather implausible. A more realistic scenario appears to be one where the star 
formation density evolves as 

M*) = 0-23 e3 , 8 f + 4 M 7 M yr -1 Mpc~ 3 . (7) 

This model fits reasonably well all measurements of the UV-continuum and 
Ha luminosity densities from the present-epoch to z = 4 after an extinction 
correction of A1500 = 1-2 mag (A 2 soo = 0.55 mag) is applied to the data pjj , 
and produce a total EBL of the right magnitude (/50 = 1). Since about half 
of the present-day stars are formed at z > 1.3 (hence their contribution to 
the EBL is redshifted away), the resulting visible mass density is fl^h 2 = 
0.0031 150 {(M/Lb)* = 8.6/50). Note that this estimate ignores the recycling 
of returned gas into new stars. 

The observed EBL therefore requires that between 7% and 16% of the nu- 
cleosynthetic baryons are today in the forms of stars, processed gas, and their 
remnants. According to the most recent census of cosmic baryons, the mass 
density in stars and their remnants observed today is fl s h — 0.00245+ , ' , l Q 2 1 g 
[ fl6| , corresponding to a mean visible mass-to-blue light ratio of {M/Lb) s = 
3.4^^3 (h = 0.5) (about 70% of this mass is found in old spheroidal popula- 
tions). While this is about a factor of 2.5 smaller than the visible mass density 
predicted by equation (]?]), efficient recycling of ejected material into new star 
formation would tend to reduce the apparent discrepancy in the budget. Al- 
ternatively, the gas returned by stars may be ejected into the intergalactic 
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medium. With an IMF-averaged yield of returned metals of y m 1.5 Z Q |] 
the predicted mean metallicity at the present epoch is yfl^j 'fib — 0.25 Zq, in 
good agreement with the values inferred from cluster abundances [p6| . 
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